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Abstract 

We obtain explicit formulas for the Neumann coefficients and associated quanti¬ 
ties that appear in the three-string vertex for type IIB string theory in a plane- 
wave background, for any value of the mass parameter /i. The derivation involves 
constructing the inverse of a certain infinite-dimensional matrix, in terms of which 
the Neumann coefficients previously had been written only implicitly. We derive 
asymptotic expansions for large p and find unexpectedly simple results, which 
are valid to all orders in 1/p. Using BMN duality, these give predictions for 
certain gauge theory quantities to all orders in the modified ’t Hooft coupling X'. 
A specific example is presented. 
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1 Introduction 


This paper continues the study of the light-cone gauge string field theory formulation of 
type IIB superstring theory in the maximally supersymmetric plane-wave geometry i-ii- 
its purpose is to derive and analyze explicit formulas for the Neumann coefficients that enter 
in the three-string interaction vertex that describes the process in which strings #1 and ^2 
join to form string #3. In particular |l], this requires an explicit formula for the inverse of a 
certain infinite-dimensional matrix called T + (y,y) (and defined below). Here, y is the mass 
parameter that appears in the plane-wave metric. It becomes physically meaningful once we 
specify the coordinate frame, since only the product P„p is invariant under a longitudinal 
boost. The combination a'P_ for the rth string is conventionally denoted ay, with the 
momentum taken to be negative for the outgoing string so that a r — 0. We make y a 
meaningful dimensionless parameter by choosing a frame for which a 3 = —1. When strings 
#1 and ^2 join to form string #3, string #1 carries momentum fraction a.\ — y and string 
7^2 carries momentum fraction (%2 = 1 — y, with 0 < y < 1. 

Previously |4j. [5|, the inverse of T + was determined in terms of a certain infinite component 
vector called Y m (y,y ) = (F^ 1 B(y)) m and a scalar function k(y,y ) = B T T^ 1 B. Here B m (y ) 
is a known infinite component vector. (The formula will be given later.) In this paper we 
obtain explicit formulas for Y m and k , and hence also for Tjl 1 . The first step is to derive a 
first-order differential equation that determines the dependence of Y m on the mass parameter 
y. Since the value of Y m for y = 0 is known from previous analysis of the flat-space problem, 
that knowledge can be used to fix the ‘initial condition’. The resulting integrated expression 
for Y m is expressed in terms of k(y, y), which still needs to be determined. Knowledge of the 
leading large y asymptotic behavior of Y m , for all m, is sufficient to completely determine 
that function. The resulting equation involves a certain integral transform. The equation 
is solved by the inverse integral transform. The particular integral transform that appears 
does not seem to be contained in the standard mathematical references (such as Bateman), 
but we have not done an exhaustive search of the literature. 

The formulas obtained at this point in the analysis are explicit and complete, but they are 
not yet in a form that is convenient for exploring large y expansions. The large y limit is of 
particular interest in light of the proposed correspondence 0 between type IIB superstring 
theory in the plane wave background and a certain sector of J\f = 4 SU(N) gauge theory, since 
the dual gauge theory is believed to be effectively perturbative in the parameter X' = 1 /y 2 . 
In the final section of this paper we present explicit formulas, which are valid to all orders in 
A' perturbation theory, but omit non-perturbative terms of order e - 27r H“d. The latter terms 
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can in principle be extracted from our results with sufficient effort. We use these formulas 
to make a specific gauge theory prediction to all orders in A'. 


2 Review of Basic Formulas 


The three-string interaction vertex for type IIB superstrings in flat space was worked out in 
H [10[ and generalized to the plane-wave geometry in |], 0. This vertex is represented as 
a state in the tensor product of three string Fock spaces, where the individual strings are 
labelled by an index r = 1,2, 3. Any particular three string coupling is then obtained by 
contracting it with three specific string states. The formula for the three-string interaction 
vertex contains a bosonic factor 


I V B ) = exp ( - ^ ^2 a rm N 


mnsn 


|o>. 


(i) 


r,s=l m,n =—oo 1=1 


The quantities jV^ n are called Neumann coefficients. The three-string vertex also contains 
a similar fermionic factor \Vp) made out of the fermionic oscillators and a ‘prefactor’ that 
is polynomial in the various oscillators. We will not discuss either of these in this paper. 
However, aside from an overall factor v(n,y) that does not involve the oscillators, they are 
constructed out of essentially the same quantities, so the formulas that will be derived here 
determine them also. 

In describing the Neumann matrices, it is convenient to consider separately the cases in 
which each of the indices m, n are either positive, negative or zero. Henceforth, the symbols 
m, n will always denote positive integers. One result of [1], for example, using matrix notation 
for the blocks with positive indices, is 


N rs = S rs - 2(C' r C'- 1 ) 1/2 A (r)T r; 1 A (s) (C' s C'- 1 ) 1/2 . (2) 


Here C mn = m5 mn and (C r ) mn = tOrm^mn are diagonal matrices, with 

UJrm = 2 + (//ay) 2 . (3) 

The definitions of the matrices and T + will be given shortly. 

The blocks with both indices negative are related in a simple way to the ones with both 
indices positive by 

= ~ (UrTT’u,) mn . (4) 

where 

U r = C~\C r - fia r ) = C(C r + /ray)” 1 . (5) 
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For the remaining Neumann 


_33 

In the case of N ' these are the only nonvanishing terms, 
coefficients the other nonvanishing terms are 


-jryrs l\f sr 
^mO = ^ Om 


= \J^- 




{C r C~ l ) l/2 A^Y^B 

a , 

r,se { 1 , 2 }, 


^oo = (—l) r+s+1 (l + yak) 


N 3r — N r3 — ~y 

iV oo ~ Jv oo ~ 1 


'OL r 


yj a r a s 

r e {1,2}, 


s,t G {l,2},r G {1,2,3}, 


( 6 ) 


where a = a.\a 2 a 3 = —y( 1 — y ), and k = k(n,y ) will be defined shortly. 

The matrices A^ r \y) and the vector B(y) do not depend on y, so they are the same as 
in flat space. Namely, 


^»=-(- ir^yfi^ yAnml,V 

7r 


n 2 — m?y 2 ’ 


(1 — ?/) sin rimy 
n 2 — m 2 ( 1 — y) 2 ’ 


7 r 

A ( 3 ) — X 
^mn — u mn, 

r _ 2 ( pm +1 sin rrnry 
m 7r y(l-y)m 3 / 2 ' 


(7) 


A number of useful identities relating these matrices are included in appendix A. Out of 
these and the diagonal matrices U r we construct 
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r + = J2 A(r)u r Air)T - 

r =1 


( 8 ) 


Note that the only y dependence enters via U r . 

The quantities that we especially would like to evaluate explicitly are the matrix Tl{ 1 (/i, y), 
the vector 

(9) 

and the scalar 

k(y,y) = B T Y- 1 B. (10) 

In the case of flat space (y = 0) the results are known. Specifically, the Neumann matrices 
in eq. (Q) may be written as 


where 


N rs = 

mn 


mna 


ma s + not. 


- N m N n for/r = 0, 


AC =- f m {-a r+1 /a r )e mTo/oir for// = 0, 


a r 


( 11 ) 


( 12 ) 
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where a 4 = aq is understood, 


fm{ 7) = 


r(m7) 


m\T(m'y + 1 — m) 


and 


r 0 = y^a r In \a r \ = ylny + (1 - y) ln(l -y). 

r= 1 

-i-l _ 1/ 


(13) 


(14) 


__ 3 

In particular, still for y — 0, = |(1 — N ), Y rn — —N m , and k = 2r 0 /a. I 11 other words, 


n,o. = 0,,)= lr f r [™ v) . e~» = 

mil (my + 1 — m) 


r(l + my)r(l + m(l-y)) ^_ mTQ 

2 T(1 + m ) 


B n 


and 


Mk = 0,!,) = —2(2^+ ln(1 ~ !/ ) 


.1 -y y 

In |Q, Q] the following identity was derived for arbitrary //, 


{r;',c 3 } = c + I 0102 


-1 , ‘ CK'YY^CK 1 . 


(15) 

(16) 

(17) 


21-)- yak 

Note that this determines T^ 1 (y,y) in terms of Y(y,y ) and k(y,y). In particular, this 
formula was shown to imply that the generalization of eq. (ED to nonzero y takes the form 


N r N s 

m n 


—rs mna 

1 -f- fiak (y. s (jj rm + OL r uj s 


where 


(18) 

(19) 


K. = - [(c- I c r ) 1 / 2 f/- 1 7 i < ’' )T r] m , 

bnt neither it nor k(y,y) was determined explicitly at nonzero y. 

Some preliminary analysis of large y asymptotics was initiated in [|2|. |3j, Q, [H|, though 
not much can be done without additional explicit formulas. It was found that the leading 
(large y) term in the expansion of TlJI 1 is given by the first term on the right-hand side of 
eq. (Uzb- Defining R by 


r ; 1 = -cc^ + r, 

it is easy to see that the leading term in R is of order y~ 4 . Specifically, 

R -> a R 7i y 2 ^~ y ^ C 3 BB t C 3 + ■■■ 
h 

where a R is a constant and the next term in the expansion is of order y ~ 6 . Similarly, 

1 dk 


k(y,y ) 


yy(l - y) n[yy(l - y)]‘ 


+ 


( 20 ) 


( 21 ) 


( 22 ) 
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Inserting these expansions into eq. (EZD, one learns that 

1 


a-R^k — 


64 ‘ 


( 23 ) 


It is very difficult to determine clr and a*, separately without additional explicit formulas. 
The asymptotic expansion of Y was found to have the structure 


Y 

1 rr, 


1 

'1 

(1 \ 

m 3 


— 

—m — 

- x 

—T • 



2 

V4 ) 

p 2 



B m . 


(24) 


The value of x is of particular interest. It was estimated numerically to be approximately 
1/16 in ||, and we will show below that this is correct. 


3 The Differential Equation 

This section describes the derivation of a differential equation involving T m (/i, y) and k(/i, y ). 
For the benefit of the reader who would like to skip the details of the derivation, and move 
on to the next section, the result is stated here: 


dY r 


ldF 


d[i [2 <9/i V VmJ wl 


Y 

l rn 


(25) 


where a; m = a; 3m = \Jm? + p 2 and 

F(/x,y) = Ml + fiak(fi,y)] = ln[l - fiy( 1 - y)k{n,y)\. 
This has the formal solution 

m ri r OF ( u \ , 

Y m (0,y). 


Ym{p,y) = -—exp 

UJm. 


1 PdFf j^ 

2 L da \ u, 1 ^ 


(26) 


(27) 


/ o M 

Thus, if we knew k(/j,,y), we would know F(y,y), and then Y m (y,y), and hence all the 
Neumann coefficients. 

The derivation of eq. (|25|) is rather involved. Let us sketch the derivation here and then 
fill in some of the details in appendix B. The matrix T + = ^ A^U r A^ T only depends on 
fj, through the dependence of U r on p. Its derivative can be written in the form 


^± = A aB B T + yN, 

Ofi 2 


where 


1 V = 


r)T 


(28) 


(29) 


r =1 
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It follows that 


(30) 


BY 1 

= -kaY - g]V-'NY. 
dg 2 


But NY can be recast in the form 


NY = 9l C^B + g 2 B, 


where the coefficients (]\ and g 2 are scalar quantities: 

2(1 + gak) 


(31) 


9i = 


2 + gak + n 2 aki 


92 = l 2 


a\ ak 2 + gakki + 2 k\ 
2 + gak + n 2 aki 


and 


h = B t C^Y. 


The above equations imply that 


dk T dY 1 2 

— = B — = -ak - gg 2 k - ngik 2 . 
o/i og 2 


(32) 

(33) 

(34) 


This is not very useful as it stands, since there is no other apparent way to determine k 2 . (/q 
could be determined, but that will turn out not to be necessary.) Substituting the equation 
for NY and an identity for [C^ 2 ,^ 1 ] deduced from eq. (|T?D, one can recast the derivative 
of Y in the form 


dY 

dg 


= (F 0 + F l C 3- 1 + F 2 Cz 2 )Y, 


(35) 


where the scalar functions F t are given by 


a 


Fn = -ak — gq 2 H— gg\ 

2 2 hy 1 + gak 

1 a 2 


(h - gk 2 ), 


F 2 — —ggi + ^k- z 9i- 


a 


-(k — gki 


Using eq. (I 


2 r 1 + gak 

for gi and g 2 , and eq. ( [34|) to eliminate k 2 in favor of k', we hnd 

a 1 


F 2 = -g, F\ = -gF 0 , F 0 = 


2 1 + gak 


(k + gk'), 


(36) 


(37) 


so that eq. fl35[) can be written in the form (|25|). 

As a first application of eq. (|25|), let us consider it for large g. The expansion in eq. ( |22| ) 
implies that at large g the leading behavior of dF/dg is — 1 /g. Substituting this into eq. (|24p, 
one deduces that x = 1/16. 
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4 An Integral Transform 


In this section we show that the known large /i behavior Y m ~ ^B rn , together with eq. (pH]), is 
sufficient to deduce an integral transform, whose solution we are able to determine explicitly, 
thereby obtaining explicit formulas for fc(/i, y), F(y,,y) and hence Y m (n,y). 

Since Y m B m for large /i, eq. (|27|) implies that 


exp <; - 


dF 


B 




B n 


2Y m (0 y 


Taking the logarithm of both sides and integrating by parts, this becomes 


(38) 


(m 2 + fi 2 ) 3/2 F(ix,y)dn = G(m,y), 


where, using eq. (|15D, 


G{z,y) = -1 o In 


r (i + z) 


z z 2 \T(1 + £t/)T(1 + z(l — y)) J 


(39) 


(40) 


Eq. (fjOp, which must hold for m = 1,2,..., and 0 < y < 1, determines F(/j,,y) and hence 
k(y,y). 

The function G(z, y) is holomorphic in the right-half z plane and goes to 0 as z —» oo 
in that half plane like a power (1 /z 2 ). These properties make this extrapolation from the 
positive integers to a continuous variable z unique. Moreover, they are exactly what is 
needed to construct an inverse transformation giving F in terms of G. Let us state this as: 
Theorem. Suppose that g(z) is holomorphic in the right-half /~-plane and vanishes like a 
power at infinity in that half plane. Then the solution of the equation 


(z 2 + x 2 ) 3 / 2 f(x)dx = g{z ) 


(41) 


(for Re(z) > 0) is given by the inverse integral transform 


f(x ) = ~i 


x 


7T 


g(—ixcosO) cos 9 d6 = —i 


x 


7T 


fw/2 


[g(—ixcos 6) — g{ixcos6)\ cos Odd. (42) 


We refer the reader to appendix C for the proof of this theorem. 
Applying the theorem to the problem at hand, we learn that 

,2 P7T 


F(y,y) = -i 


■ B 


71 


G(—iy cos 9, y) cos 9 d6. 


Substituting the formula for G, 


2i 


F(y, y) = 2 ft 0 H— 


r dd ln 

T(1 — iy cos 9) 

/o cos 9 

T(1 — ifiy cos 6*)T(1 — ifi{ 1 — y) cos 6) 


(43) 


(44) 
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The formula for Y m also requires the derivative 


— -^ — - = 2r 0 H— f dd[0(l — iy cos 0) — ?/0(l — iyy cos 0) 




-(1 - 2 /) 0 ( 1 - *A*(1 - y) cos 0 )], 

(45) 

where, as usual, 0 = (lnT)'. 




Substituting the expansion 




0(1 + Z) '■ 

X ^ 

- 'y +z '£ n{z+n) 

n=l 

(46) 

and using the integral 

s>7 r 

dO 7T 

(47) 


,/o a + 6cosd y/a?-b 2 

gives the result 

dF(y,y) 

dy 

oo 3 

n =1 r =1 UJrn 

(48) 


Since F(0,y) = 0, this integrates to 

F(y,y) = 2r 0 /i + 2 ^ In [(cu rn + /xa r )/n] = 2r 0/ u - 2 lndet(f/if/ 2 f/ 3 ). (49) 

n, r 

Note that det(t/i?72^73) is convergent even though the individual det(17 r ) diverge. We can 
regulate them in a way that does not change the product (since ^ a r = 0) by recasting F 
in the form 

F(y,y) = 2^11 + 2^2 (f)(ya r ), (50) 


where 


00*0 = ln 


/ n 2 + x 2 + x \ x 


oo 

£[■ 


arc sinh 


x\ x 


n/ n 


The regulated determinant of U r is then exp[— <j>(ya r )\. Substituting the expansion of dF/dy 
into the formula for Y m gives 


Y m (y, y) = exp (y - u m )r 0 + J^(0 r - 0 mr ) 

r=l 

where 0 r = (p(ya r ) and 


ln 


“ 1 “ 



where, we remind the reader, uj m = \Jm 2 + y 2 and u rm = \Jm 2 + (ya r ) 2 . 

The formulas © and (|52|) are explicit expansions of F and Y m which converge for all 
finite y. Thus, in a sense, they solve our problem, ffowever, they are not yet in the most 
convenient form for exploring large y expansions. 


5 Asymptotic Expansions for Large // 


In this section we develop asymptotic large y expansions for all quantities that appear in 
the three-string vertex. It turns out that the quantities all have essential singularities at 
y — oo arising from terms proportional to e - 27r /4«r-|. The existence of such terms, which 
correspond to non-perturbative effects in the dual gauge theory, has been noted in ||. We 
proceed under the assumption that y\a r \ is sufficiently large (for all r) that these terms 
can be neglected, and we use the symbol ~ to denote this approximation. We stress that 
the formulas we present encapsulate all orders in a power series expansion around A' = 0 
from the dual gauge theory point of view. However, they do not match smoothly with the 
known flat space results at /r = 0 because the omitted terms become important in this limit. 
Finally, the limit y —> 0 (at fixed y) is also of interest, because it can be used to extract 
vertex operators for the emission of on-shell particles [|9|]. Our asymptotic formulas are not 
well-suited for studying this limit because we omit terms of order e~ 2nfiy . 

Let us demonstrate how to use the integral equation (|39|) directly to fold the constant a k 
introduced in eq. (|22|) . It follows from eq. (^2|) that 


F(y, y) - - In [nfiy(l - y)/a k } + • • ■ (54) 

at large y. Let us define F(y, y) = F(y , y) + In [7ryy(l — y)ja k ]. Then eq. (|39|) implies that 

'mny(l-y)' 


(m 2 + y 2 ) 3 ^ 2 F(/i, y)dy = G(m, y) H-- In 






(55) 


Now let m = riX and y = A w. After scaling out A, we fold 


(w 2 + n 2 ) 3 ^ 2 F(\w,y)dw = A 2 G(\n,y) + In f^ 7m y^ 


A 2 n 2 


2 Clk 


(56) 


In the limit A —> oo the left-hand side goes to zero, while from eq. (|4(J) it is easy to check 
using Stirling’s approximation that the right-hand side goes to ^ ln(4afc). This determines 
a k = 1/4 and, from eq. (|23|) , = A. 
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5.1 F(fjL,y) and k(fj,,y) 


Let us proceed by studying the function cp(x) defined in eq. (fTTp. Since </> is clearly odd, it 
is sufficient to consider large positive x here. Taking two derivatives gives 


(x = 




1 1 
+ 


' (x 2 + n 2 ) 3 / 2 x 2x 2 ’ 


(57) 


where we have used eq. (|146|) to evaluate the sum (the exact result is given in eq. (|156|) 
Integrating twice with respect to x leads to 


4>(x) ~ — (x + |) lnx + C\X + c 2 , 

where c* are constants of integration. Inserting eq. (|58|) into eq. ( |50| ) gives 

F(v, y)~~ Hwi 1 - y )] + 2c 2 . 


(58) 


(59) 


The constant c\ has dropped out since ^a r = 0, and we can determine c 2 = hi 47T by 
comparing eq. (|59|) to eq. (^4|) with — |, obtaining thereby 


F(l*,y) ~ - \n[Anfiy(l-y)]. 


We emphasize that eq. (p0|) is much stronger than eq. 


(60) 


when we wrote the latter, we 


might have expected corrections involving powers of l//i, but in eq. ((6^) we have proven 
that the only corrections are exponentially small. An alternative derivation of this result is 
presented in appendix D, where we show that the absence of power law corrections follows 
from a simple contour integral argument. In appendix F we derive the exact formula 


F(fi,y) = — ln[47r/«/(l — y)] + J(fiy) + 


where 


J(x) = 


ln(l 


0 —2 7TXZ\ 


7T J l Z\Jz 2 - 1 


-dz. 


(61) 


(62) 


It is easy to read off all the exponential corrections to eq. (|6l]) by writing out the series 
expansion of the logarithm. However, we will not keep track of these exponential corrections 
in the following sections, and instead simply use eq. (|60|) and the definition (|26l) to write 

1 


k(y,y) 


l 


yy0--y) ^^ 2 y 2 (i-y) 2 ' 


(63) 
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5.2 Y m (n,y ) 


Although it is straightforward to develop an asymptotic expansion for the function </>„ 
dehned in eq. (|53|), a more direct route is simply to rewrite eq. (|7]) as 

"1 f°° OF (_ ii \ , dF 


m 


Y m (y,y) = —exp 


UJr] 


dy 


1 - — 1 dy - i — 
0J m J ./„ oy 


1-— 1 dy 

UJ, 


Y m (0,y). (64) 


The first integral is just eq. (|38|), and the second integral is elementary after substituting the 
asymptotic expansion F' « —1 /y. The hnal result is 


Y m (y,y) 


y + uj m m 1 \j y + + h 2 

J^m. — 1 / 777 /, 


2/r 2u; rn 2\A2y -^/m 2 + y 2 

which can be conveniently summarized as 

Y 


1 


_=[/ 3 1/2 C' 3/2 Cj' 1 i?. 

2v/2^ 3 3 


(65) 


( 66 ) 


5.3 The Neumann vectors 

The final step in the construction of the Neumann matrices involves evaluating the matrix 
products A^ t Y which appear in the Neumann vectors in eq. (jl9i). From eq. (^5|) we have 


{A^ T Y) n 

Using eq. we find 

(A (r)T Y') f 


- 1 ) 


77-Li /- OO 

^ \/n v—^ sin 


E 


in 2 (7rmy) \Zh + V d 2 + m2 


y n 2 y 2 (l - y ) ^ mn 2 - n 2 /y 2 yj^ + m 2 


:(-l re {1,2}. 


ny( 1 - y)^V^y 

The Neumann vectors in eq. (|l^) can therefore be expressed as 

1 


(67) 


( 68 ) 


N' 


n: 


2vr?/(l - y) 


{-l) r{n+1) y/o^(2ymj rn Urn) 1/2 , re {1,2}, 


— — (2yU3 n U3n) X ^ 2 B r 


Actually we can combine these expressions in a useful way. If we define 

Sim = s 2 m = 1, s 3m = -2 sm(nmy ), 

then we have simply 

_„ ^_Xj r U+ 1 ) 


n: 


\Oi r 


2ny(l -y)\ 2 ynuj rn U r . 


re {1,2, 3}. 


(69) 


(70) 


(71) 
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5.4 Consistency checks 


We should check that the vector Y we have found indeed satisfies B T Y « k and r + y ps B. 
First we have 


B t Y 


1 oo 

1 v—sm 


7T 2 y 2 (l -y) 2 2^2jl^ i m 


E 


2 ( 7 r my) y /' + C/' J — 






(72) 


With the help of eq. (|151|) , it immediately follows that this expression equals eq. 

In 0 it was shown that r + = 2C 3 C~ 1 — H, where H is given up to exponential corrections 
by 


H, 


mn o o 

/i z 7T z 


(— 1 )" l+n y/mn sin(7 rmy) sm(nny ) / dz 


V* 


1 (z 2 + m 2 / fi 2 )(z 2 + n 2 / jj?) 


(73) 


(The integral is easily evaluated, but it is convenient to leave eq. ( |73| ) in this form for the 
calculation.) The condition T + Y = B which we would now like to check is equivalent to 


HY = 2C 3 C~ l Y - B. 
Using eqs. (0) and (|65| ) we can write 

2V2n 2 B n f 


(74) 


(HY)r 


7 T 


Vz 2 - 1 ^ sin(7T my ) 2 \J d + W 1 + 


rn 


■y/Jl J 1 z 2 + n 2 /ii 2 ^m 2 + z 2 n 2 y/ 777,2 + ^2 


(75) 


After substituting eq. ( |152| ) for the sum, the remaining integral over z is elementary and 
takes the form 


Re 


dz 


V — 1 + iV z 2 — 1 

z(z 2 + a 2 ) 


7T 


y/2 c 


1 - -^Ji + ViTV 2 

V 2 v 


(76) 


for a = n//n Assembling all factors from eqs. ([75]), ([A]) and (|152|) , we find 

tin -),,« -d„ i--Ly/i + Vi + »V^ 


(77) 


Recalling eq. (|65|), we see that the desired relation eq. ( |74| ) is indeed satisfied up to expo¬ 
nential corrections. 


5.5 Some remaining quantities 

The Neumann matrices are completely determined, to all orders in l//r, by the factorization 
identity (|T8D and the Neumann vectors given in eq. m ■ For the sake of completeness, we 
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catalog here two more quantities of interest, which follow easily from our results. The first 
is the matrix R defined in eq. (|20|). By comparing the expressions (|2j) and ([T8|) for r = s = 3 


we can determine T + 1 , and hence R, in terms of Y m . Using eq. (|66|), we arrive at 


Rmn = (T 


- 1 \ 


m 


2 vr 


mn(-l) 


m+n 


sin(7 rmy) sin(7 my) 

^n)'s/U rn U n 


(78) 


where (as before) 

Uj m = sjm 2 + p 2 , U m = — + - . (79) 

m 

Finally, we remarked above that the prefactor is polynomial in quantities A" and Y which 
are constructed out of bosonic and fermionic oscillators respectively (see |2], ||, 0]). The 
normalization of X and Y involves a factor called /(/i) in [Sj]. It may be obtained from the 
formula |2], [J 

/(At) = \/-27r« 2 «3 lim (-^(Cd^r^B)^ (80) 

n—>oo 

It follows immediately from eq. (|B5|) that 


/ 0 ) 


y/4irfj.y(l-y) 


\/l + fiak. 


(81) 


In fact, closure of the supersymmetry algebra requires f(fi) = + fiak 0. Note that this 

/(/i) is separate from a still undetermined function u(/i), to appear below, which appears as 
an overall factor in the cubic part of the Hamiltonian and supersymmetry generators. 


5.6 Summary of Neumann matrices 

We summarize here the final expressions for the Neumann matrices obtained in this paper. 
As before, we use the notation 

ai=y, <*2 = 1 -y, a3 - -1. ( 82 ) 


Sim = s 2m = 1, s 3m = -2sm(nmy), 
and u rm = \Jm 2 + (/ray) 2 . Then for m, n > 0 we have 


(83) 


-^frs _ 1 (-l) r ( m+1 )+ s ( n+1 ) j\a r a s \[u} rm + fia r ){uj sn + fia s ) ^ 

™ mn ~ ~ A / Srm.Sfj 


27T OL s UJ rrn “j - Ot r UJ sn 


(jJrmfjJ. 


(84) 


rmy-' sn 


—rs 1 (-i)dm+i)+d«+i) ha r a s \{u rm ^ fia r )(u sn - fia s ) 

2 V _ m _ n r ''-' _ \ / S rrn S .qn • 


27 T Q. s UJ rrn “I - Oi r LxJ sn 




(85) 


rm^sn 
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As before, the symbol m denotes that we have omitted terms of order e 27 I A a d (for all r). 
For r = s = 3 these are the only nonzero components. If we define 


1 

s io — S 20 — y=> 


S30 — 0 , 


( 86 ) 


then eq. (|84|) continues to hold when m is zero and n is positive (or vice versa). Finally, if 
n — rn — (1 we have 



1 (-l) r+s 
A-rji y /a r a s ’ 


r,s 6 {1,2}, 


(87) 


_Q r _ r q _ 

N a> = N m = re {1,2}. (88) 

We remark here that we have been assuming throughout this work that /i is positive, 
since it is clear from the analysis of [![] that only the absolute value of /i enters into the 
bosonic matrix elements. The behavior of various quantities under /i —» — /i was exploited in 
Q to derive several useful identities. Although we have not considered this involution here, 
it may be worthwhile to do so. 


6 A Matrix Element 


In this section we use our result to calculate, to all orders in A', a particular matrix element 
of the Hamiltonian which has so far only been computed to first order in X' in the dual gauge 
theory i m 


. We then explain in detail how this matrix element is encoded in the gauge 


theory. 


6.1 String field theory 

Consider for m, n > 0 the three states^ 

( 1 | = |(01 + 

< 2 | = < 0 |, 

(31 = |(01«-ifo_ n )K+ m J _ n ), (89) 

where i and j are SO(4) indices. The two-impurity states (1| and (3| decompose into the 

1, 6 and 9 representations of SO(4). For definiteness, we fix ? ^ j, and we could choose to 

symmetrize or antisymmetrize in ij at the end of the calculation. Actually, it turns out that 

the matrix element vanishes when either two-impurity state is in the 6. 

1 We caution the reader that the basis of oscillators employed here and in 00 differs from that used by 

H and most gauge theory papers by the transformation a® MN = ^j(o|„| — i sign(n)a_| n |) for ti/ 0. 
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We have not discussed the prefactor in this paper, but it has been shown in pi] that 
for states of the form (^) (in particular, for states with no fermionic excitations), the three- 
string coupling | H) in the Hamiltonian is given effectively by 


\H) = v{»,y)V\V B ), 
where | Vb) was defined in eq. (Jl|) and we have defined^ 


OL \ ^ 


r=1 


E ^rm / \ f 

- e{m)al m a r 

llOL r 


m= 1 


e(m) = 


m > 0 , 
m < 0 , 


(90) 


(91) 


where the SO(4) index (which we have suppressed) is contracted between a' and a. Finally, 
the function v(y,y) is a measure factor which has not yet been determined. It arises from 
the path integral which defines the cubic string vertex. In flat space, the function can 
be determined by Lorentz invariance (it is 1 in the supersymmetric theory, and some very 
complicated function of y in the bosonic theory). The plane wave superalgebra does not 
have enough generators to fix this overall factor, although comparison with gauge theory 
requires that v(/j,,y) —> 1 for large y. The corresponding factor in the supergravity vertex in 


the plane wave background has been determined to be a constant |ll 
For the states (^) we find the matrix element 


Hnmy = (l|(2|(3|tf) = v{fJl, y) 


1 a 
4 2 


^lm ^3 r 


7 13 


2 iV mn -2 lV_ m ,_„) 


7 13 


Using eqs. 


and 


y,ot 3 

, we hnd after dramatic cancellation the simple result 
sin 2 (7r ny) \v(ii,y)/j> 2 y 


H ss -(1 

nmy 2 /r 2 1 


y)- 


7 T- 


(92) 


(93) 


^lm^3n 

The quantity in brackets is equal to 1 at leading order for large /i, reproducing the result 
presented in [BJ. As is now standard in the literature, it should be understood that the 
cubic interaction H considered here enters the full Hamiltonian with a coefficient equal to 
the effective string coupling, g 2 = ingsy 2 . 


6.2 Relation to gauge theory 


The result (93) provides a concrete all-loop prediction for the gauge theory, which we now 
explain. Consider (in the large J limit) the normalized BMN [| 8 ], [14|, 15, U| operators 

1 J 

O n = V p 2vikn / J TrU l Z k M Z J ~ k ), 


2 The apparent discrepancy between this formula and the one given in Q is entirely due to the change of 
basis in footnote 1. 
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Jy 


TV = 1 1 

" y/Jy( 1 - y ) ^JW+z ^ 


^ e 2nikny/j . Tr(0 < zV^ Jl ' _fc )Tr(Z (1 ~ I ' )J ):, 


yy _ 




( 94 ) 


v/jv 7 ^ 

where i / j G {1,2, 3,4} are S0(4) indices (see the discussion below eq. (j39|)) labelling four 
of the six scalar fields of A/” = 4 SU(N) gauge theory, and Z = ^(</> 5 + i</> 6 ). 

At zero string coupling (g 2 = 0), the single- (double-) trace operators O (T) defined 
in eq. (|94|) correspond respectively to one- (two-) string states and have definite conformal 
dimensions 


A n — J + 2y /1 + X' 


n* 


A y n = J + 2^1 +\'n 2 /y 2 , 


A y = J + 2. 


(95) 


For finite g 2 , these operators mix ]1J, [15], [l(|, and the state-operator correspondence has 
been worked out to order g 2 in || [H| (see also [0). The operators which correspond to the 
desired one- and two-string states are 


O p = O p 


92 

2 ^ 


k= — C 50 ' 

oo 


dy C pky T» - | / dy C„T « + O(gl), 


where 


T t = T k ~ -j °pky°P + ^(triple trace) + 0(g$), 

p=— OO 

1 ~ V sin 2 (7r py) sin 2 (7r py) 

Jy 7r 2 (p — k/y) 2 ’ py ^JiT 2 p 2 ' 


C p ky 


(96) 


(97) 


The triple-trace operators in eq. (|96|) will not be important for this calculation. The anoma¬ 
lous dimension matrix element between the single-string state O and the two-string state T 
is read off from the two-point function 


(2Trx) An+Aym (O n (x)Tl(0)} = -g 2 h mny ln((rA) 2 . 

The prediction from (|93[), when expressed in gauge theory variables, is simply 


hmny 


H, 


mny 


_ A 11 

y/Jy^-v)™ 2 


y sin" ( 7 my) 


[1 + An 2 ]" 1/2 [l + \'m 2 /y 2 ] 


/ 2 /„ ,21 - 1 / 2 


Jy 


7T 


(98) 


(99) 


up to the overall function v(fi,y) discussed in the previous subsection. The leading X' term 
in this result agrees with the one-loop field theory calculations of |T3], 0 m when the 
appropriate operator redefinition in eq. (^6;) is taken into account. Note that since eq. ( |99| ) 
is already proportional to A', probing the subleading terms in this expression would require 
a two-loop gauge theory calculation, which has not yet been reported in the literature. 
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7 Conclusion and Discussion 


The primary results of this paper are twofold. First, we have shown through a quite intricate 
analysis that it is possible to determine all Neumann coefficients in the plane wave back¬ 
ground exactly in terms of a single function F(fi,y), for any value of the mass parameter /i, 
and we have provided an explicit formula for F in eq. Secondly, we have investigated 

the large /i behavior of these coefficients and presented simple formulas in section 5.6 which 
give the Neumann matrices to all orders in a l//i expansion. Note that although we have 
not discussed the matrix elements of the prefactor (nor the fermionic Neumann matrices) 
in detail, they are very easily obtained from the Neumann vectors (|7T|) using the results of 

a a an- 


Perhaps the most remarkable fact about our results is that although the Neumann ma¬ 
trices are very complicated functions of fi, the final expressions (|84|) and (^) for the 1 /y 
expansion are very simple. Ours was a long (and perhaps circuitous) road, and although we 
have derived a number of nontrivial identities along the way, one cannot help but wonder 
whether there is a more direct path which yields the same final result. In particular, our 
prediction ( PU| ) for a particular gauge theory calculation is so simple that it cries out for 
explanation by some clever argument, perhaps along the lines of [[IT]. 

The origin of this simplicity is the fact that the function F(y,y) behaves for large fi like 


F(y,y) = -In[ 47 r/iy(l-y)] + 0(e" 2 ^, e" 2 ^, e“ 2 ^M), (100) 


with no perturbative corrections (i.e., inverse powers of /i). This remarkable fact, which we 
have proven in two different ways (in appendices C and D), is reminiscent of various non¬ 
renormalization theorems. Indeed, although the BMN operators in eq. fl94|) are non-BPS in 
general, their anomalous dimensions 

A(O n ) - J - 2 = 2Vl + A'n 2 - 2, A' = ( 101 ) 

J z 

are nevertheless finite in the limit of large’t Hooft coupling, provided that J is simultaneously 
taken to infinity so that X' remains finite. This suggests that there should be some residual 
‘effective supersymmetry’ protecting these operators and their interactions. 

We have not worked out explicit formulas for the non-perturbative terms in eq. ( |100|) , 
although these could in principle be obtained by extending the analysis of the appendices. 
Finally, the presence of fractional powers of A' in certain string field theory observables 
has been noted in 0 0 These surprising powers appear in matrix elements when the 
total number of ‘impurities’ is not conserved. For example, a matrix element in which two 
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impurities are created or destroyed is in general a factor of fi = 1 /x/A larger than a similar 
matrix element with conserved impurity number. Although this seems enigmatic from the 
dual gauge theory point of view, it has been stressed |3], [L 8 | that there is no reason why the 
large J limit of the A expansion has to agree, order by order in A, with the X' expansion - 
especially in light of the fact that the BMN limit requires taking A —> oo, which need not 


be a trivial extrapolation. For the anomalous dimensions eq. (' ll Oil) it seems to work, but in 
order to reproduce impurity non-conserving interactions on the gauge theory side one might 
have to sum the A expansion to all orders and then take the large J limit. 
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A Basic Identities 

The infinite matrices Amu, C and the infinite vector B rn satisfy a number of useful relations 
which we record here: 


A( r)T C'A (s ) = — C5 rs , r, s G {1,2}, 

(102) 

a r 


A {r)T= a r C~ l 5 rs } r, s G {1,2}, 

(103) 

A^ t CB = 0, rG {1,2}, 

(104) 

b t cb= . 2 , 

y( 1 - y ) 

(105) 


where we have used <23 = —1, Some additional useful identities are 

3 1 

V — A (r) CA (r)T = 0 , (106) 

Ot T 
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( 107 ) 


^cv4 M C _1 .4 (r)T = - BB r , 

r =1 - 

where, as before, 

a = a 1 a 2 a 3 = ■— 2/(1 — y)- (108) 


B Derivation of the Differential Equation 


In this appendix we fill in the details of the derivation in section 3, making frequent use of 
the identities presented in appendix A. 

We start by deriving eq. fl3ip . Using eq. (|102|) to simplify T+CA^ gives the identity 


Y- l C :i A (r) = CA (r) - ^-T-'A^Cr, r G {1,2}. 


It follows from this that 


Y T CU 3 1 A^ = -—Y T A {r) U; 1 C ) r G {1, 2}. 


a 


(109) 


( 110 ) 


Multiplying eq. m on the right by A hi and using eqs. (|109|) and (|110|) then gives 


T^A^Cr - OrCsTl'A^ = -- 

+ 3 + 21 +flak 


CU 3 1 YY T A^U~ 1 C, re {1,2}. (Ill) 


Comparing this to eq. (0). we can write the following formula for all r, 


r l 1 A ir) C r - a r C7 3 r-M (r) = C5 r3 - cu^yy^a^cu- 1 . 

2 1 + fiak 


( 112 ) 


Multiplying this further by B T C 3 1 on the left and by a r C l C r 1 fih) T on the right, then 
summing over r, gives 


-ahB T - Y t N = -B T C 3 2 - -- 

2 3 2 1 + fiak 


(.k — /j,ki)(-akB T + fiY T N). 


(113) 


The transpose of eq. (|113|) is a linear equation for NY whose solution is ([31]). 

Next we obtain eq. fl35[). We start by using a variant of eq. ( |T7| ) to obtain an identity for 
[c 3 2 , r; 1 ] = which yields 


1 ry 

[<+, r;']B = K fc - Mi)(<+ - M+) - (h - Mn)(i - MV 1 )]*" (114) 


and hence 


r; 1 c 3 - 2 fi = 


a 


2 1 + fiak 


{kofi — k\) + (k — fi z k 2 )C 3 1 + (—b kfi + k\f+)C 3 2 


a 


Y. (115) 


Then we can substitute eq. m into eq. (|30|) and use eq. (|115|) to arrive at eq. 
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C On an Integral Transform 

In this appendix we show that the general solution of the Fredholm integral equation of the 
first kind 



(116) 


IS 



[g{ix cos 9) — g(—ix cos 9 )] cos 9cl9. 


(117) 


For the amusement of the reader we present two derivations of this result. The first is a set 
of manipulations that indicate how we arrived at the solution, while the second is a more 
rigorous proof of the theorem as stated in the text. 

C.l An elementary manipulation 

As integral equations of the first kind are notoriously difficult to solve, let us attempt to 
circumvent the problem. The method used here does not seem to be in the canonical texts 
and is hoped to be of some use. 

Consider the following integral identity p0|: 




(118) 


Next suppose that zg(z) can be expanded as a power series in e z — 1: 



(119) 


Therefore by eq. (|118|), we have 




from which we can directly read off the solution 



( 120 ) 
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(a) 




(b) 


(c) 


Figure 1: The contour C for appendices C and D. In (a) and (b) we show the singularity 
structure of eq. (|123|) , with branch points at —1 and —y'/y. In (c) we show the structure of 
eqs. ( |126| ) (for integer k ) and ( |140|) . 


This formidable sum can actually be performed due to the identity 


from which we obtain 


J 0 {kx) = - / e tkxcos0 d9, 

7T ./n 


oo ] 


j =0 k =0 

oo 


S(x) = 


/ 7r 

J2aj(e icosxe -l) j d9 
j =o 


— 


7T 


2 /"IT 


d9 cos 9 g(—ix cos 9). 


( 121 ) 


( 122 ) 


Therefore all dependence on the coefficients cij drops out and we have a generalized Fourier 
transform of rather simple form. 

The main shortcoming of this derivation is that it assumes the existence of a rather 
peculiar expansion, which might not be necessary. The alternative approach presented in 
the remainder of this appendix is clearer in this regard. 


C.2 Representation of a delta-function 

Here we consider the integral 

x(y,y’) = j— I -/== 7 — f)3/2 , ( 123 ) 

4?u J c a/1 + w [wy + y') s ' 2 

where y,y' > 0 and C is the contour shown in fig. 1. For y' < y the singularity structure is 
shown in fig. 1(a), and the contour can be pushed off to infinity, giving zero for the integral. 
For y' > y, the contour encloses no singularity (fig. 1(b)), so the result is again zero. 

Since x(y,y') vanishes for y' ^ y , let us check whether it is a delta-function (and in 
particular, that it does not involve any derivatives of delta-functions) by integrating it against 
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the test function e ty ': 


dw 


/ dy' x(y,y') e ty — 7T~ / /- __ 

/o 27T? J c a/ 1 + w-y/u; 


[l + e twy ntwy(<$>(\J twy) — 1)] , (124) 


where the result of the y' integral involves the error function 


2 r 

*(*) = -= / 
do 


; 1 dt. 


Next we make use of the elementary integral 

1 f dw 


-_w k = 


(125) 


(126) 


27t i J c a/1 + w^/w v k 
In particular, note that this is zero when k is a positive half-integer since the branch cut then 


runs from —1 to — oo, so there is no singularity within the contour. The term in eq. (|124|) 
proportional to “—1” therefore integrates to zero, since it contains only half-integer powers 
of w. For the remaining terms we use the expansion 


00 O/u r^k 

Zj 


1 + e x y/nx^(y/x) = 7 ^ 7 — 


k =0 


(2k — 1)!! 


(127) 


Note that (—1)!! = 1. Combining eqs. ( |127| ) and (|126| ) leads to 

°° 2 k (ty) k 


/»00 ^ 

/ dy' X (y,y')e~ ty ' = Y 
d° k =o 


(2k — 1)!! V k 


1 \ oo 

2 I _ 


E ^ 


k =0 


k\ 


(128) 


where we have used the identity (2k — 1)!! = (—2 ) k k\( 2 ). If x contained any terms propor¬ 
tional to derivatives of delta-functions, we would have obtained a polynomial in t times e~ ty 


in eq. (|128|). Since this did not happen, we have proven that 


x(y,y') = 


\[y 


dw 


= d(y-y'), y,y'> 0. 


(129) 


47T* J c yfl+w (wy + y’) 3 / 2 

More precisely, we have proved that / x(y,y') f(y')dy' = f(y) for any / that can be written 
as a convergent Laplace transform. It is not excluded that an even weaker condition would 
suffice. 


C.3 Formal proof 


To verify our solution let us substitute f(x) from eq. (|117|) into the integral equation (|116|) 
and replace the integration variable x by t — xcos 6. This gives 


i 

71 


r /2 t 2 [g (-it) - g (it)} cose _ _i_ 
J 0 (t 2 + z 2 cos 2 9) 3 / 2 n 


dt 


t[g(~it ) - g(it)] 
(t 2 + z 2 ) 


(130) 
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where we have used 


r*7r/2 


cos 9d9 


J 0 {f + cos 2 Of/ 2 t(t 2 + z 2 )' (131) 

We can now recast the remaining integral in eq. (|130|) as a contour integral with the contour 
C enclosing the positive u axis. This gives 

g(-iy/u) 


- / d u ' 

2 m J c , [u + z z ) 


= s(z), 


(132) 


where we have used the assumed analytic and asymptotic properties of g, which are just 
what are needed to deform the contour into one that encircles the pole at — z 2 , thereby 
reducing the integral to exactly g(z). 

We have proven that the function f(x) given in eq. (EH, when substituted into eq. ( |116|) , 
indeed gives g(z). It remains to show that the result is unique, in other words that when 
g(z) from eq. (EH) is substituted into eq. (|117|), we recover f(x). To this end, we define first 


f{x ) = xf(x 2 ), 

and let v' = x 2 so that the integral equation (EH) becomes 

- / dv'(z 2 + v')~ 3l2 f{v') = g{z). 


(133) 


(134) 


Our proposed solution is, in terms of f(v), 


f(v) = 


i\ v 


7T 


g(—i\/v cos 9) cos 9d9. 


Letting yTc = —i cos 9 gives 




(135) 


(136) 


where the contour C is as in fig. 1. Substituting eq. ( fL 16| ) for g(z) and using the delta-function 
representation eq. ( |129| ) establishes eq. ( |136| ) and hence completes the proof. 


D Asymptotic Behavior of F(/i 1 y ) 

The analysis of appendix C was completely general, and we only used rather weak assump¬ 
tions about the form of g(z). Here we use the contour integral techniques above to study 
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the particular case of interest, with g(z) given by G(z,y) from eq. (]40j). Recall Stirling’s 
asymptotic expansion 


lnr(l + z) w (z + |) \nz — z + - In 27T + ^ 


y 2n+l ' 


(137) 


n =0 


The coefficients c n may be expressed in terms of Bernoulli numbers, but they will turn out 


to be irrelevant. The crucial fact is that only odd inverse powers of z appear in eq. ([137]). It 
follows that 

1 .22, h (iA 

(138) 


G(z,y) =3 “ln(2irzs(l - y)) + 

n= 1 


2n+l ' 


Now letting y/w = —icos 6 in (|43|), as in the previous subsection, gives 

■ 2 r dw 


F(y,y) = 


y~ 


-.G(yy/w,y). 


(139) 


(140) 


27 t 7 J c y/l + W 

Consider first the inverse odd powers in eq. ( |138| ). They lead to integrals of the form 

r dw 1 
Jc \/l + w W n+l / 2 

for n a positive integer. The branch cut is now as shown in fig. (lc), so the contour may be 
pulled off to infinity, giving zero. Therefore only the first term in eq. (|138|) contributes, so 

1 f dw 


F(y,y) 


2 ni J c w\/l + w 

\w{2-Kyy{l- y))F— 


\w{2-k y\/wy{l - y)) 


dw 


In w 


zdw. 


(141) 


2 m Jc w\/l + w 47 ri J c w\/l + w 

The first integral gives 1 since it just picks off the pole at w = 0. It is an interesting though 
straightforward exercise to check that the second term gives — In 2. This completes the direct 
proof of eq. 


E Advanced Sums 

Let us review an elementary trick which can be used to evaluate certain infinite sums. 
Consider first the sum 


/M’ /( 


z = 


Z 2 + X 2 V 2 y/z 2 + X 2 ’ 


for x > 0 and v > 1. This sum can be written as the contour integral 

1 


2 ni 


dz f(z) 7T cot 7 rz, 


(142) 


(143) 
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Figure 2: This figure shows the contours relevant for analyzing the various sums and integrals 
in appendix E. In each case, the contour is deformed from the one enclosing the real axis 
to the one enclosing the branch cuts on the imaginary axis. Note that the integrands in 
eqs. (|144|) and (|152|) have two additional poles on top of the branch cuts, while the integrand 


in eq. (|148|) has two additional poles on the real axis. 


where C passes from — oo — ie to +oo — ie slightly below the real axis, and then returns 
slightly above the real axis (see figure 2). The contour may be deformed away from the poles 
of f(z)7r cotTTZ on the real axis to pick up the other singularities instead. In this case the 
only other singularities are branch cuts from ±ix to ±foo, and poles at z — ±ixv sitting on 
the branch cuts. Rescaling z and combining the two cuts, we find therefore that 


£ 


n =—oo 


l l 

n 2 + x 2 v 2 sjn 2 + x 2 



dz coth(7ra;z) 
Vz 2 - 1 £ 2 - v 2 


(144) 


where P stands for the principal value. This is an exact formula, but if we are not concerned 
with terms vanishing exponentially for large x, then we can set coth(7nc,z) = 1 in the integral, 
which is then easily evaluated to give 


£ 


1 1 

n 2 + x 2 v 2 y/ n 2 + x 2 


2 arc cosh(u) 
x 2 v \Jv 2 — 1 


(145) 


As in section 5, the symbol ~ denotes that we have dropped terms of order e 2nx . Taking 
the limit v —> 1 from above finally gives the result 


£ 


1 

(■n 2 + x 2 ) 3 / 2 


2 


x z 


(146) 
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For sums in which f(z) itself contains trigonometric functions, it is convenient to expand 
these and write the sum in the form 

OO 

£ [/,(«)+ (-l)”/ 2 (n)]. (147) 

n =—oo 


The term with fi(z) can be evaluated using the above trick, while the second term can be 
evaluated by deforming the contour integral of f 2 (z)7T esc irz. Applying this method and 
taking into account both the contribution from the poles at Am/y and the discontinuity 
across the branch cuts gives 


£ 

m =1 


sin 2 (7 vmy) \Jx + \fx 1 + m 2 
m 2 - n 2 /y 2 sj x 2 + m 2 


1 „ f°° dz y/l + i\fz r ^\ 1 

2x 3 / 2 J 1 — 1 - n 2 / (x 2 y 2 ) F(x, y, z) 

(148) 


when 0<?/<l,x>0, and n is a non-zero integer. We have defined the function 


F(x,y,z) 


1 

2 


[coth(7r xyz) + coth(7rx(l — y)z)\. 


(149) 


Ignoring exponentially small corrections, we can set F — 1 and evaluate the integral in 
eq. (|148p explicitly, arriving at 


OO 


sin 


! ( nmy ) \Jx + \Jx 2 + m 2 ny 3 1 2 \j \/ n2 + x2 y 2 "%V 


m 


m 2 - n 2 /y 2 p x 2 + 


m z 


An ybi 2 + x 2 y 2 


(150) 


The case n — 0 must be considered separately since the poles at m 2 = Azn/y are then lost. 
The result in this case is 


OO 

v—> sm 


in 2 (7r my) \Jx + \fx 2 + m 2 7t 2 ?/(1 — y) 


7T 


m =1 




Vr + 


m 


V2x Ay/2x 3 / 2 ' 


(151) 


A combination of all of the above techniques is needed to tackle the final sum 


£ 

m=l 


sin 2 (7r my) \Jx + \[x 2 + m 2 
m 2 + x 2 v 2 \/ x 2 + m 2 


JL RpP r dz Vl + I 

W 2 ./, VFWI 2 2 -t< 2 

7T Im v^l — i\/u 2 — 1 

’ 


(152) 


which is valid for u > 1. 
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F Some Exponential Corrections 


In this appendix we derive the exact formula (0) which allows systematic determination of 
the exponential corrections to F(n,y). It follows from eq. (|144|) that 



n =1 


(153) 


where 


I(x) = P J 


00 COth(7 TXZ) 


(.z 2 - l ) 3 4 5 / 2 


1 

2 


: 



Integrating the second term by parts gives a divergent piece that cancels the divergent piece 
from the first term, leaving 



(155) 


Therefore 



00 zdz 1 


(156) 


V z 2 — 1 sinh 2 (7 txz) 


Integrating twice with respect to x as in subsection 5.1 and using eq. (|50|) yields the formu¬ 


las (|6l|) and (|6^). 
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